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I. INTRODUCTION 



Persistent equilibrium current occurring in isolated mesoscopic normal metal rings pen- 
etrated by an Aharonov-Bohm (AB) flux is one of the most interesting phenomena in meso- 
scopic systems It is periodic in the flux with period of a flux quantum (or half a flux 
quantum) and has either diamagnetic or paramagnetic sign according to different experi- 
mental situations 0,0]. It has been widely accepted that spin-orbit (SO) interaction may 
play a crucial role in the orbital magnetism, e.g. the sign and the period of currents. Levy et 
al. have suggested that persistent currents observed by them may have a diamagnetic sign 
due to the effect of strong SO interaction (although they pointed out that their determina- 
tion contains some ambiguity) 0. On the other hand, Altshuler et al. have claimed from 
thermodynamical arguments that the current averaged over spatial disorder is always para- 
magnetic even in the strong SO interaction limit Q. Extensive studies done subsequently 
T0| , [TTH have confirmed the result of Altshuler et al., and furthermore have revealed some in- 



teresting aspects for SO effects on the orbital magnetism. Particularly in a non-perturbative 



approach employed in [TT| universal and nonuniversal aspects of SO effects have been dis- 
cussed, such as the reduction factor of currents, etc. The above interesting studies on SO 
effects, however, have been concerned with a free electron model without electron-electron 
interaction. It is hence desirable to examine how electron-electron interaction is combined 
with SO interaction to affect persistent currents in an AB geometry. 

In this paper we wish to investigate the effect of SO interaction on persistent currents 
in mesoscopic rings of mutually interacting electrons for the canonical ensemble. For this 
purpose, we study the Hubbard ring with SO interaction by combining the non-perturbative 
treatment of SO interaction of Meir et al. [12| with the Bethe ansatz technique [|l^,|l^]. We 
find that the interplay of SO interaction and electron-electron interaction produce some new 
effects on the orbital magnetism. We further point out that simple reduction factors of 
currents due to the SO effect are modified in the presence of electron-electron interaction. 

The organization of the paper is as follows. In Sec. II, we briefly depict how to diagonalize 
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the Hubbard ring with SO interaction by the Bethe ansatz method, and derive the excitation 
spectrum exactly. The key point is that the non-perturbative treatment of SO interaction 



for a non-interacting case |12| is still applicable to the Hubbard model because of local 
SU(2) symmetry of the Hubbard interaction. We then study, in Sec. Ill, the SO effects on 
persistent currents for case of the canonical ensemble, putting stress on the sign and the 



period. As is the case without SO interaction |T3|,|i"5|], the current shows a quite different 



behavior according to the number of electrons N c modulo 4. We will see that upon averaging 
over strong SO interaction, some new effects on persistent currents are brought about by 
the electron-electron interaction. In Sec. IV, a possible extension to more general models 
including long-range interaction is presented. Summary and conclusion are given in Sec. V. 

II. HUBBARD RING WITH SPIN-ORBIT INTERACTION 

We consider the mesoscopic Hubbard ring with SO interaction. The Hamiltonian reads 

H = ]T UiS^cl^c^ + h.c. + c h c n4i c n- (!) 

1,(7,(7 L 

where {Si) aa > is a matrix of SU(2), and ti is chosen to be real for the sake of time reversal 
symmetry. The spin-dependent hopping matrix in ([!]) reflects the effect of SO interaction. 
The interplay between the Hubbard interaction and the SO interaction makes it difficult to 
treat the model directly. It is found, however, that one can still deal with the Hamiltonian 
exactly by combining the Bethe-ansatz technique with the non-perturbative approach de- 
veloped for a non-interacting case [12|. The point is that the hopping term in (jl]) with SO 



interaction is cast into the diagonal form in spin space by a unitary transformation 

= jj(i) Sl . . . S ._ u ( 2 ) 

with the matrix being chosen appropriately. The transformation rotates the frame of 
spin space by a different angle at each site so as to produce the homogeneous hopping matrix. 
An important point is that the onsite Coulomb term in ([!]) is invariant under such a local 
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spin rotation because of its SU(2) symmetry |16|]. Hence one can formally gauge away SO 
interaction, and find the conventional Hubbard model in the new spin space, 

i,a i 

where a denotes a transformed spin variable, and f and j label "up" and "down" spins re- 
spectively in the new spin frame defined differently on each site. The effect of SO interaction 
is now incorporated into spin-dependent twisted boundary conditions for the eigenfunctions 



m 



fe{N + 1) = exp[ ? ($ + &Wff(l), (4) 

where exp(i&5) is the eigen value of Si • ■ • Sjy, a = ±1 corresponds to "up" and "down" 
spins in the new spin space, and $ is an AB flux in unit of $o/(27r) with the flux quantum 
$ — h/e. The electron wave function acquires an additional phase shift ±5 due to SO 
interaction after transversing the ring. We will refer to S as the SO phase shift hereinafter. 
It is noteworthy that the above trick to simplify the model is not specific to the Hubbard- 
type interaction, but is also applicable to any interaction, including long-range type, which 
has local SU(2) symmetry. Such extensions will be discussed later in this paper. 

In the following, t n is assumed to be site-independent so that the lattice is regular without 
spatial disorder. Thus we put £j = —t. It is now straightforward to obtain the Bethe 
ansatz solution to the Hamiltonian (JJ) with twisted boundary conditions @ |T7 |. Following 



a standard method, two kinds of rapidities are necessarily introduced to diagonalize the 
Hamiltonian. For the ring system with N sites, one thus gets the coupled transcendental 
equations for charge (pj) and spin (A a ) rapidities JT^JTj 



N s 

Pj N = 2TxIj + ($ + 5) + J2 vi^Pj ~ Xp), (5) 

13=1 



X V {sm Pj - X a ) = 2tt J a - 25 - £ v((^ - A a )/2), (6) 

3=1 f3=l 
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with i](p) = — 2tan~ 1 (4tp/£7), where the number of total electrons (down-spin electrons) is 
N c (N s ). The quantum numbers Ij and J a are integers (or half integers), which specify 
charge and spin excitations. The total energy is given in terms of the charge rapidity, 
E = —2tJ2f=i cos Pj- Note that the effect of SO interaction effectively shifts the quantum 
numbers Ij and J a whereas that of Coulomb interaction appears via the two-body phase 
shift function t](k). 

Let us now consider the effect of SO interaction on the energy spectrum. Applying a 
machinery developed by Woynarovich |13[ to Eqs.([5|) and (P), one can readily classify the 



excitation spectrum including the finite-size corrections which are important for the meso- 
scopic Hubbard ring. For the fixed number of electrons, low-lying excitations are specified 
by two kinds of quantum numbers D c and D s , which respectively carry the momentum 
4kpD c (charge current) and 2kpD s (spin current), where kp is the Fermi momentum. The 



excitation spectrum including SO interaction is written down as [|13| , |19 



where v c and v s are the velocities of charge and spin excitations respectively, and K p is 



the critical exponent for the Akp oscillating piece of charge correlation functions [13]. For 



the Hubbard model 1/2 < K p < 1 [[H5|-|2"T|] . These three fundamental quantities charac- 



terize the Luttinger liquid properties of interacting electrons completely P^] , which can be 
straightforwardly evaluated using the Bethe-ansatz integral equation resulting from (5) and 
(6) |T3|JT9| -^T| . For a given number of electrons it is necessary to find the lowest energy 
state correctly in order to derive the expression for currents. It is crucial for this purpose to 
notice that the quantum numbers D c and D s respect the following selection rule reflecting 
the Fermi statistics [|l"3|j , 



N c + N s + 1 

D c = ^ (mod 1), (8) 



We note that in the absence of SO interaction, the above spectrum has been analyzed by 
Yu and Fowler to study persistent currents p4 ]. 



III. EFFECTS OF SPIN-ORBIT INTERACTION ON PERSISTENT CURRENTS 

We now study the effects of SO interaction on persistent equilibrium currents in the 
Hubbard ring. Since we wish to discuss the magnetism for the canonical ensemble with the 
fixed number of electrons, let us briefly summarize some results known for the canonical 
ensemble. The property of the orbital magnetism is sensitive to the number of fermions 



carrying the currents p3| , |24| . For example, a free electron model without SO interaction the 
orbital magnetism depends on the total number of electrons N c modulo 4; i.e., the ground 
state is diamagnetic for N c = An + 2, paramagnetic for N c = An with period of a flux 



quantum, and paramagnetic with period of half a flux quantum for N c = An + 1, An + 3 [15 



The results are slightly modified in the presence of electron-electron interaction as shown for 



the ID Hubbard model []14| , Pq| . For instance the paramagnetic state for N c = An is altered 
to a diamagnetic one by electron-electron interaction except near half filling (one electron 
per lattice site). In the presence of SO interaction, further modifications are expected to 
occur in the orbital magnetism due to the interplay of SO interaction and electron-electron 
interaction. In particular we will see below that for a certain parameter regime of the 
interaction strength, the paramagnetic state is stabilized for N c = An + 2 in contrast to the 
case without SO interaction for which the ground state is always diamagnetic. 

In order to clearly see what is going on, we first study the case of An + 2 in detail, 
and mention the other cases later in this section. It may be plausible to introduce here 
an important key quantity v c K p /v s , which will be helpful for following arguments. For 
noninteracting electrons, v c K p /v s = 1 for any electron concentrations because v c = v s and 
K p = 1, whereas in correlated cases of U ^ the value of v c K p /v s ranges from to oo 
depending on the interaction strength as well as electron concentrations. We show v c K p /v s 
as a function of electron concentrations in Fig. 1. Roughly speaking, one can see from 
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this quantity whether the SO effect may be enhanced or suppressed by the electron-electron 
interaction. For example, the SO effect is suppressed for v c K p /v s > 1, whereas it is enhanced 
for v c K p /v s < 1. 



A. Persistent currents for N c = An + 2 (N s = 2n + 1) 



It is seen from Eqs. (^j) and (|9|) that the selection rule for the quantum numbers in this case 
is D c = (mod 1) and D s = (mod 1), which implies that the ground state in the absence 



of SO interaction is diamagnetic around $ = PD| , |23|P^| . Turning on an AB flux, there 
occur two different situations according to the magnitude of v c K p /v s . It is straightforward 
to derive the persistent current / = —dE/d§ from Eqs. (|7|). For v c K p /v s > 1, the current 
takes the form |14] 



ttN v 



7T 



rrN 



where 



7T 

$c= 77 + 



-7T < $ < — $ c 
-$ c < $ < $ c 
< $ < 7T 



(10) 



2 ' 2K p v c 



K p v c 



(11) 



One can see that the SO phase shift 6 simply alters the critical value <3> c of the AB flux, 
and the diamagnetic nature around $ = is not modified by SO interaction. From Eq.(O) 
it is explicitly seen that the effect of SO interaction is suppressed as the value of v c K p /v s 
increases. Such cases with v c K p /v s > 1 realize at lower electron concentrations for U ^ 0, 
as seen from Fig. 1. It is instructive to point out here that for 5 = 7r/2 the effects of 
electron-electron interaction disapear in Eq . (|TT|) , and hence the period is halved as has been 
known for a free electron model fT2|. An alternative expression of Eq . (|10"D in Fourier series 
expansion is found to be more convenient for following discussions, 

-l) n 



/($) = £ v ' K " 



iiN n 



cos 



71V 



IT 

n\ — 

2 2K p v c K p v c 



, v s 5 



sin(n$). 



(12) 



In contrast to the above case, the expression for currents in the case of v c K p /v s < 1 
depends on the value of the SO phase shift 5. For < 5 < nK p v c /2v s + tt/2 the current is 
given by the same expression as Eq. ([370]) , whereas for ttK p v c /2v s + 7r/2<5<7ritis changed 
to 



-^($ + 7T) -7T<$<0 

ttN 

-^($-7T) 0<$<7T. 

ttN 



(13) 



This is cast into an alternative formula in the Fourier series expansion as 

W=E^sin(n<&). (14) 

One can see from Eq . (|13"D that the current has a paramagnetic sign around $ = 0. This 
change in sign is due to the interplay of the SO interaction and electron-electron interaction, 
which still plays a crucial role upon averaging over the strong SO interaction (see below). 



B. Strong SO interaction limit 

We have seen that the SO phase shift 5 crucially modifies the characteristic behavior 
of the orbital magnetism in correlated electrons. The quantity 5 is to be determined in 
the range [0, ir] for a given scattering process by SO interaction. For example when the 
scattering length of SO interaction would be comparable to the length of the ring, the value 
of 5 could be around tt/2, which results in the period of currents halved by the effect of 
the SO interaction as seen from (|T2|). In more general cases when the length of the ring 
may be much larger than the scattering length, one has to take into account all possible 
rotations in spin space due to SO scattering. Such a limit is referred to as the strong SO 



interaction limit pl| , p!2| . Since the SO phase shift 5 is directly related to the rotation angle 
in spin space, one should average currents over 6 from to n with the weight sin 2 5 |TTLP7 . 



The Fourier transformed formulae of currents are more convenient to carry out the average. 
Using Eqs. flT2|) and (|T4"D, we thus derive the expression for persistent currents in the strong 
SO interaction limit. In the case of v c K p /v s > 1 (low electron densities), we get 
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It is seen that all the odd harmonics are dropped, and hence the period is reduced to half a 
flux quantum. In contrast to a free electron case in which only harmonics n = 0, 2 remain 
and the harmonics n = 2 gives a diamagnetic sign change with a reduction factor 1/2 ||11|| , 
all higher even harmonics survive in the correlated case. However it is found that the ground 
state is still diamagnetic around $ = 0. To see this more explicitly we evaluate the sum of 
the coefficient by Fourier transformation: 

^f v c K p v s v s K p m \ f v s m7i \ _ 1 

^{mv s + 4v 2 K 2 -Av 2 m 2 ) [ L) Sm [v c K p J~ 2' ^ 



which results in the diamagnetic current = — (v c K p /7rN)$>. Note that this is exactly 
same as the current without SO interaction. Hence we can conclude that in case of v c K p /v s > 
1, the behavior of currents around $ = is not modified by the SO effects even in the strong 
interaction limit. This is because the SO interaction can affect the occupation of energy levels 
only when there exists a finite AB flux in the system. 

The situations are somewhat different for v c K p /v s < 1 (close to half filling). There are 
both diamagnetic and paramagnetic contributions to currents depending on 5. The current 
average over the strong SO interaction is thus given by 



n 
+7T 



vtn* . fTxv r K n \ %v^K^ . / nix nirv 9 



v r .K n v 2 n 2 . / i , 

sin H 1 — 7, sin — + 



2n 2 Nn Y^v 2 c K 2 - v 2 n 2 \ v s ) v s n(Av 2 K 2 - v 2 n 2 ) \ 2 2v c K f 
irv c K, 



x sin(n$). (17) 



The period is not halved upon averaging over 5 in this case. For small positive $ with 
1/3 < v c K p /v s < 1, one can perform the inverse Fourier transformation to get for $ — > 0, 



v c K p /tcv c K p \ v s -v c K p 
cos 



v a \ v. 



2 



TlV c K p . f7TV c Kp 

7T — sm - 

V. V Vs 



(18) 



which has a paramagnetic sign. Since paramagnetic contributions to the current averaged 
over S increase compared with diamagnetic ones when v c K p / v s is decreased, we can say that 
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the current for < v c K p /v s < 1 always shows a paramagnetic sign around $ = 0. We would 
like to stress here that this paramagnetic state realizes as a consequence of the interplay of 
SO interaction and electron-electron interaction, and should not appear if either of the two 
is absent. 

We have computed persistent currents using the formulae (^3j) and flT7|) together with 
the Bethe equations for v c ,v s , and K p p|. In Figs. 2 and 3, persistent currents in the strong 



SO limit are shown as a function of the AB flux $ for two different cases. 



C. Reduction factors in currents 



Here we make a brief comment on reduction factors in persistent currents by the SO 



effects. According to Meir et al. |T2[ and Entin-Wohlman et al. [[LI], reduction factors due 
to the SO interaction are written in a simple and universal form for a noninteracting model 
with even number of electrons. An essence of their idea is that the persistent current 
in the presence of SO interaction can be expressed in the form 

= !„($ + 8) + /„($ - 8) (19) 

where Iq is the current per each spin component without SO interaction. This leads to 



universal reduction factors in the form of Fourier expansion [11 .13], 



/($) = ^cos(n<5)a n sin(n$) (20) 

n 

with a n being the Fourier coefficients for the case without SO interaction. One can see that 
the reduction factor cos(n<5) depends only on the SO phase shift. This expression leads to 
a rather simple result that upon averaging in the strong SO limit, the current is reduced 
by a factor of 1/2, and the fluctuations by a factor 1/4 [p], 12|. For interacting electrons, 
however, it is seen that the SO phase shift should not appear in a simple form of $ ± 8. It 
is instead combined with a factor of v c K p /v s reflecting the electron correlation effects. Thus 
in the presence of electron-electron interaction, the reduction factors follow from Eqs.(|12"D, 
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COS 



7T 

n[ — 



TtV. 



2 2v c K p 



+ 



v c K p 



I cos 



7T 

n\ — 



7TV a 



2 2v c K p 



(21) 



for v c K p /v s > 1. For v c K p /v s < 1 we can not define the reduction factor in such a simple 
form, as seen from Egs- fllTf ). So, the above remarkable properties obtained for the SO 
effects based on a free electron model are changed when the electron-electron interaction 
is introduced. To avoid confusions we would like to mention that the present nonuniversal 
results may not be contradicted with the universal reduction factor of currents expected 
in disordered system [ rU| , |IT , p? -30|. In such cases the average over disorder may play an 
essential role, which has not been taken into account in the present calculation. 



D. Electron- number dependence 

We have been concerned so far with the case of N c = An + 2. As mentioned before 
the results are sensitively dependent on the number of electrons. Here we summarize the 
results for other cases. The calculation can be performed in parallel to the above example 
of N c = An + 2. 

(a) N c = An (N s = 2n). In this case the selection rule for the quantum numbers reads, 
D c — 1/2 (mod 1) and D s = (mod 1). For v c K p /v s > 1 the persistent current in the 
Fourier transformed form is then given by 

1 ; ^TiNn I \2 2v c K p v c K p )\ y h y ' 

which shows a diamagnetic sign around $ = 0. In the case of v c K p /v s < 1, where only the 
paramagnetic state realizes without SO interaction, the expression is the same as Eq.(p2|) 
for 7r/2 — v c K p 7i/2v s < 5 < it, while Eq.(|T4|) for 0<<5<7r/2 — v c K p 7i/2v s . Thus the para- 
magnetic state near the metal-insulator transition (half-filling) is changed to a diamagnetic 
one by SO interaction for tt/2 — v c K p 7r/2v s < 5 < tt. It is remarkable, however, that in the 
strong SO interaction limit we obtain exactly the same results in currents as for the case of 
N c = An + 2 upon averaging over all possible phases 5 (see Figs. 2 and 3). 
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We can hence summarize the results for the even number of electrons as follows; the 
ground state is always diamagnetic for v c K p /v s > 1 (low electron densities) in the strong SO 
limit, whereas for v c K p /v s < 1 (near the metal-insulator transition point) the paramagnetic 
state is stabilized not only for N c = An but also for N c = An + 2 as a result of the SO effects. 

(b) Odd case (N c = An + 1, An + 3). We first consider the case for N c = An + 3. Since 
Eq.([?p holds only for the case with iVj > iVj_, the roles of "up" spin and "down" spin for 
$ > are interchanged for $ < in Eq. (]?]). Thus the result for $ < in the case of odd 
number of electrons can be deduced from that for $ > by changing 5 — > —5. Noting the 
selection rule D c = 1/2 (mod 1) and D s = 1/2 (mod 1) for $ > 0, the current in the ground 
state is now given by 

V c K p ( _ 7T 



,,,$ + - -7T < $ < 

vriVV 2) - (23) 

v c K p f^ ^ 



$ < $ < 



7T 



ttN V 2, 

which has a period of half a flux quantum, and exhibits a paramagnetic sign around $ = 0. 
It should be noted that the above expression is independent of 5, and hence there are not any 
modifications due to SO interaction. The effect of SO interaction indeed appears in the next 
order corrections in 1/N, which are not taken into account in our formulation. For a non- 
interacting case the next-order corrections have been evaluated by Entin-Wohlman et at, 
resulting in a small shift of energy minima from $ = ±7r/2 to $ = ±[(1 + l/N c )ir/2 — S/N c ] 
n|,|3T|, which may be actually neglected in mesoscopic rings. We note that for another odd 



case of N c = An + 1 , we get the exactly same results of currents as for N c = An + 3. 

In summary we can say that characteristic properties of persistent currents in the strong 
SO limit are classified by the parity of the electron number even for correlated electron 
systems. 



IV. EXTENSION TO MORE GENERAL MODELS 



In the Hubbard model the electron-electron interaction is assumed to be short-ranged 
(on-site), which in turn enables us to treat the model exactly. In more general cases, however, 
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effective interaction would be of long-range type, since the screening effect of the Coulomb 
interaction may become less effective for mesoscopic metallic rings. In such cases it is quite 
difficult to get the energy spectrum exactly. So it is desirable to find a possible way to extend 
our analysis to more general cases including long-range interaction. We wish to briefly depict 
a simple idea how to treat such cases. 

We recall here a trick used to simplify the Hubbard model with SO interaction, i.e. a 
unitary transformation which incorporates the SO interaction into the boundary effects @. 
Note that this technique is still applicable to more general long-range interactions so long 
as they retain local SU(2) symmetry. Such a local SU(2) symmetry for interaction may be 
expected to hold in ordinary cases, such as partially screened Coulomb interaction, etc. The 
remaining problem is then how to obtain the expression for the low-energy spectrum like 
Eq.([?D including the SO effect and the AB flux. To this end the bosonization technique 
may be useful [^,^,^,^] because it can formally describe low-energy states even for non- 
integrable systems. The bosonization scheme has been previously used by Loss to discuss 



persistent currents for a spinless fermion system [24 



Following a standard technique in bosonization p2| , |33| , we now discuss how the SO effect 
on currents can be treated in nonintegrable systems. We do not have to specify an explicit 
form of interaction here, but only assume the interaction V(r) to be invariant under local 
SU(2) transformations. In general, low-energy gapless states of ID metallic electron systems 
compose of two independent Luttinger liquids ]22| corresponding to charge and spin degrees 
of freedom [|T9H2l|| . Hence the system is described by the sum of two Gaussian models with 
conformal charge c = 1. 

Let us now introduce the boson fields for spin and charge degrees of freedom. It is 
found that the spin-dependent twisted boundary conditions (|2|) due to an AB flux and SO 
interaction are incorporated into boson phase fields as p2]j24] 



7T 2 j. TCI* 

<t>o = £ TFZ e [°L + + <Po, P + M PT , (24) 



13 



E 



7i 



2Lk 



e ikx [a\ a + a^ a ] + + M a 



TCX 



E 



7T 



^ = E 



2Lk 
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sgn(fc)e 



ifor f t 



k,p 



2Lk 



'' sgn{k)e lkx [al a - a_ v ] + V + + v 7 ^) I + f ) 



(25) 



(26) 



(27) 



where ak )P (a k p ) and afc )fl -(a][ (7 .) are boson annihilation (creation) operators for charge and 
spin densities, respectively, and M p (M a ) and J p (J a ) are the charge (spin) number and 
the charge (spin) current, respectively, which are defined by M p = (Mf + Mj_)/\/2, = 
(M T - M{)ly/2, J p = (J T + Ji)/V2, and J ff = (J T - JJ/V2. We note here that M T(i) and 
J"m) describe topological excitations introduced by Haldane p2|j . They satisfy the selection 
rules, (— l)^ ^ 1 = (_i) M s+ J s ; w h ere s — j ; | anc i iVofn) is the total number of up (down) 
spins. 4>o,u and 9 ^ are conjugate variables of J v and M u , respectively. 

Using the above phase fields we can write down the low-energy effective Hamiltonian, 
from which the persistent current directly follows, as long as the system belongs to the 
universality class of Luttinger liquids. The energy spectrum of the effective Hamiltonian 
with finite-size correction terms reads [f22 



71 

2L 



v pM M 2 p + v pJ J p + V2- + v uM Ml + v oJ [Ja + V2- 



(28) 



TT J \ 71, 

where v u m and v u j (y = p, a) are Luttinger liquid parameters describing the velocity of 
excitations. It is to be noted that the effect of electron-electron interaction V(r) is only to 
renormalize these parameters. The ground state is given by the condition, M p = M a = 0, i.e. 
Mj = Mi = 0. Thus using the selection rules mentioned above we obtain the topological 
constraints for Jj and Jf (1) Jj even, J^ even for N c = 4n + 2, (2) J^ odd, J^ odd for 
N c = 4n, (3) Jj odd, Jj even (and vice versa) for N c = An + 3, (4) J\ even, Jj odd (and vice 
versa) for 4n+ 1. Consequently, if we put v p j = K p v c and v a j — v s , we find that the energy 
(P8|) with these topological constraints is equivalent to Eq.(^), and D c and D s correspond 
to Jf/2 and (Jj_ — J\)/2, respectively. 
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Using Eq.(p8|) and the above topological constraints, we can obtain the persistent current, 
and discuss its sign and period quite similarly to the case for the Hubbard ring. Therefore 
all the expressions derived in the previous section can be directly applied to the present 
case by regarding v c K p /v s as a free parameter to be determined. In order to determine 
the parameter, we generally need input data from another microscopic calculations, e.g. 
numerical diagonalization. We note that several elegant techniques to get the Luttinger 
liquid parameters numerically for nonintegrable systems have been already developed |22],|20 
which will be helpful for us to evaluate the persistent currents explicitly. 



V. SUMMARY 

In this paper we have discussed the effects of SO interaction on persistent currents in the 
mesoscopic Hubbard ring. We have investigated the problem exactly combining the Bethe- 
ansatz solution with a unitary transformation which incorporates the SO effects into spin- 
dependent twisted boundary conditions. It has been shown that characteristic properties of 
the orbital magnetism in the Hubbard ring are classified according to the value of v c K p /v s 
and the number of electrons modulo 4. In particular, we have demonstrated that v c K p /v s is 
an important key quantity to see whether the SO effects are enhanced or suppressed by the 
electron-electron interaction. 

In the strong SO interaction limit it has been found that the formula obtained for currents 
is classified by the parity of the electron number. For the even number of electrons, the 
ground state is diamagnetic with period of half a flux quantum for v c K p /v s > 1 (low electron 
densities), and paramagnetic with period of a flux quantum for v c K p /v s < 1 (close to half 
filling). In particular the paramagnetic state for v c K p /v s < 1 is realized by a combined 
effect arising from the interplay of SO interaction and electron-electron interaction. In the 
Hubbard model, the condition v c K p /v s < 1 is satisfied near half filling which implies that 
the system would be close to the Mott insulator. Therefore such a novel phenomenon for 
v c Kp/ v s < 1 is expected in general to occur for interacting electrons in a metallic phase close 
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to the Mott insulator. In contrast to the even case, the ground state for the odd number 
of electrons is found to be always paramagnetic with period of half a flux quantum, which 
is not affected by SO interaction as long as the corrections up to the order of 0(1 /N) are 
concerned. 

In conclusion the effect of SO interaction together with electron-electron interaction gives 
rise to a novel and qualitative change in the orbital magnetism for ID interacting electron 
systems. There remain several important problems to be investigated. For example we 
have not considered the effect of disorder in this paper, which would induce interesting 
phenomena together with the SO effects as well as with the correlation effects. Furthermore 
an extension of the theory to multichannel cases is desirable to confront the results with 
various experiments. These problems are now under consideration. 
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FIGURES 



FIG. 1. 

Plots of v c K p /v s as a function of electron densities n for the Hubbard model. The half 
filling corresponds to n — 1. 

FIG. 2. 

Persistent currents plotted against $/27r for U/t = 4 and n = 0.65 (v c K p /v s = 1.32) in 
the case of even number of electrons. The current normalized by t/N is shown. 

FIG. 3. 

Persistent currents plotted against $/27r for U/t — 4 and n = 0.95 (v c K p /v s = 0.58) in 
the case of even number of electrons. The current normalized by t/N is shown. 
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